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Abstract 

We study the Schwinger mechanism for the pair production of fermions in the presence of an 
arbitrary time-dependent background electric field E(t) by directly evaluating the path integral. We 
obtain an exact non-perturbative result for the probability of fermion-antifermion pair production 
per unit time per unit volume per unit transverse momentum (of the fermion or antifermion) from 
the arbitrary time dependent electric field E(t) via Schwinger mechanism. We find that the exact 
non-perturbative result is independent of all the time derivatives d , where n = 1,2, ....oo. This 
result has the same functional dependence on E as the Schwinger's constant electric field E result 
with the replacement: E — > E(t). 
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Using a proper time method Schwinger derived the following one-loop non-perturbative 
formula 

dW e 2 E 2 ~ 1 

= > — e i eB i (1) 

d 4 x 4vr 3 ^ n 2 1 J 

for the probability of e + e~ pair production per unit time per unit volume from a constant 
electric field E via vacuum polarization [1]. The distribution of the e + (or e~) production, 
d 4 xcp PT ' can no ^ °e obtained by using proper time method and a WKB approximate method 
was used for this purpose in 2j. A path integral method can also be used to obtain the pt 
distribution |3|. In the case of fermions in QED one finds 
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dW \eE\ T t 

Log[l-e *na*r]. (2) 



d i xd 2 pT 47T 3 

Recently this approach has been extended to QCD in SU(3) gauge group by directly 
evaluating the path integral and the Pt distribution of the quark and gluon production 
rate are found to depend on two independent Casimir invariants: C\ = [E a E a ] and Co = 
[d a b c E a E b E c } 2 where E a is the constant chromo-electric field with color index a=l,2,..8 3j. 

In a recent paper Q] we studied, for the first time, the Schwinger mechanism in the pres- 
ence of arbitrary time dependent electric field E(t) by using a newly derived shift theorem 
{5I for the case of scalar field theory In this paper we will extend this study to fermions. 

We obtain the following non-perturbative formula for the probability of fermion- 
antifermion pair production per unit time per unit volume per unit transverse momentum 
(of the fermion or antifermion) production from arbitrary time dependent electric field E(t): 

dW \eE(t)\ T _ T 4+- 2 , 

-Log[l - e i efi ( f )i ]. (3) 



d A xd?pT 47T 3 

This result has the remarkable feature that no time derivatives d (where n=l,2,...oo) 
occur in the final expression. This result has the same functional dependence on E as the 
Schwinger's constant electric field E result with the replacement: E — > E(t). 

Now we present a derivation of eq. (jSJ). 

The Dirac Lagrangian density for fermion in the presence of a background electro- 
magnetic field A^(t, x, y, z) is given by 

C = - e4) - m]V> = $M[A\ ip (4) 
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where = The vacuum to vacuum transition amplitude in the presence of background 

field A is given by 

The one loop effective action becomes 

S**- 1 - 1 = — «Trln[(^ — e^f) — m] + iTrln[^ — m\. (6) 
Since the full trace Tr is invariant under transposition we also get 

S-W = -iTrln[(^ - eJjC) + m] + iTrln$ + m] . (7) 
Adding eqs. ([H]) and ([7]) we find 

25 (1) = -iTrln[(p - eA) 2 + -a^F^ - m 2 } + zTrln[p 2 - m 2 } (8) 
where the Dirac tensor <j fiv is given in terms of Dirac 7 matrices as follows: 

^ = ^[7Y-7Y]- (9) 

Since it is convenient to work with the trace of the exponential we use 

a = r°° ds [e „ ls(b „ l€) _ e -is(a-ie)l (10) 
JO s 

Hence eq. ([HD becomes 

2^(1) - ^Tr — [ e - is liP-eA) 2 +!la^F liv -m 2 ~ie] _ e -is[p 2 ~m 2 ~ie]^ ^-Q 

The trace Tr is given by 

Tr = tr^troirac (12) 

where 

tr^O = f d 4 x < x\0\x > . (13) 



We assume that the time dependent electric field E(t) is along the z-axis. We choose the 
Axial gauge A 3 = so that only 

A = -E(t)z (14) 
3 



is non- vanishing. Using eqs. (jl~2]) . (fi~3l) and (fi~4")) in eq. (TTTT) we obtain 

f + OO 



2 jo s 
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dt < t\ dx < x\ 

oo J —oo 



is[(po+eE(t)z) 2 — p 2 — •pi r — m 2 +i'y 'y i E(t)— ie] — is(p 2 — m 2 — ie)l 



dy < y\ dz < z\ 

J —oo 

z > \y > \x > \t > . 



(15) 



Inserting complete set of \px > states / d 2 px\pT >< Pt\ = 1 we find (we use the normal- 



ization < q\p >= ^= e tqp ) 



SO) 



5 Tr Dirac [°° - [ d 2 x T I d 2 p T e^T +m2+ ^i[ +00 dt < t\ f °° dz < z\ 

JO S J J J —oo J -co 



2(2tt) 



-is[(-i^+eE(t)z) 2 -p 2 z +n°^E{t)] 



z > \t> - dt dz 
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(16) 



Unlike the constant electric field calculation ji], 6], the Dirac trace in the above equation 
can not be performed outside the "£" trace because the Dirac matrices are now multiplied 
by the time dependent electric field E(t). The Dirac matrix 7] 



7°7 3 



a 



(T 3 

a 3 



10 
-1 
10 
0-100 



is symmetric and hence it can be diagonalized by an orthogonal matrix U: 



(17) 



(TV) = f/[(7°7 3 )Dia g onal]f/- 1 . 



(18) 



To perform the Dirac trace in eq. (1161) we diagonalize the 7 u 7 d matrix and find: 



(7 7 ) Diagonal 



Hence we find from eq. (1161) 



Ax 

A 2 

A 3 

A 4 



10 
10 
0-10 
-1 



= v r — 1 d 2 x T i d 2 p T e is ^T +m2+ ^[( +o ° dt < t\ r°° 

2(271") ~{Jo S J J J-oo J-oo 



(19) 



dz < z\ 



Z > \t> 



dt / dz- 



Atts 



(20) 
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At this stage we use the shift theorem 



dy / dxfr(y) < x\ e - My)x+ ^ +b{ ^ )+c{y)] \x > f 2 {y) 

J — oo 

= / dy rdxfM < x - -LA | e-^ 2 ^^ \x - > fi{y) 

J J-oo a{y) dy a(y) dy 

(21) 

to obtain: 



s m 



<t\ f + ^dz<z + ^^\e-^ e2E2 ^ 2 -^ +lX ^\z + ^^ >\t>- fdtfdz ' 



oo J — oo 



E(t)dr E(t)dt 1 y 7 4vrs J 

(22) 

where the z integration must be performed from — oo to +oo for the shift theorem to be 

n 

applicable [5j. 

Inserting complete set of \p z > states (using / dp z \p z >< p z \ — 1) as appropriate we find 

4 f°c ds 



2(2tt) 



2 E/f 7 / d 2 p T e ls{p T +m2+le) 



/+oo /■ /■ r+oo ^ 

.oo * < j| y y dq * l dz k z + W)di lp > > 

< >< *|* + 4^ > |« > - /* /* (23) 

Using < >= -±== e tzPz we find 



2- 

1 



c(i) = ; 

" 2(2tt) 2 



£j jf°° L | rf 2 XT | rf 2 p T e is{ ^ +m2+ " ) [F,(s) - Jdtjdz^-} (24) 



where 



OO J J J —oo 



Fj(s) = -j^-r I dt<t\ I dp z I dq z I dz e ^e-^)* Pz < Pz | e «[-e a ^W* a 4^-iA i s(t)]| & > 



1 d 



e W)-t q * e - izqz \t > . (25) 

1 d , 

It can be seen that the exponential e B W d * contains the derivative 4 which operates on 
< p^e' 81-62 ^^^"'^'^^ > e~W)^ qz hence we can not move e ~"^)^ Pz to right. 

We insert complete sets of \z > states and \p Q >states using the completeness relations: 

/oo roo 
dz\z >< z\; 1= dp \po >< Po\ (26) 
-oo J —oo 
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to find 

r+OO 



F ' [s) = W) 



dt J dpo J dp' J dp ! Q J dp" J dp z J dq z J dz\ J dz 2 J dz 



< t\p >< Po \e izp *e-^* Pz \p' >< p' \ < p t \ Zl >< Zl \ e ^- 2 E 2 (t)z 2 +Pi-i^E(t)]^ 2 >< ^ > y, > 

< pZ\e^£ q *e- izq *\p™ >< p"\t > 

l r+oo r r r r r r r r r+<x> 
= j—^ J dt J dp J dp' J dpi J dp" J dp z J dq z J dz 1 J dz 2 J dz 

jtp 0e izp z < po | e -^y>,| p / > e ~iz lPz < p / Q | < ^isl-e^E^zHpi-i^Eit)]^ > | p // > 

e iZ2qz < pl\e^^ qz > e - izq *e- itp ° . (27) 

The matrix element, < zi\e ls ^" e2E2 ^ z2+ '^'" lXjE ^ \z% > can be evaluated in terms of the 
eigenstates of the Harmonic Oscillator with time dependent frequency. The Hamiltonian 

F=^ + ^W] (28) 

can be diagonalized at every time t in terms of time dependent creation and annihilation 
operators 



m = ,m z+ « Pt , „t (t) = ,m z - « P ,. (2 9 ) 



'2u(t) V l \/2u(t) 

Hence we find 

H = u(t)[b\t)b(t)+ 1 -]. (30) 
The eigenvalues of the Hamiltonian are just 

Xn = u(t)[n+-], n = 0,1,2,- •-. (31) 

and the normalized eigenstates < z\n t > are the usual harmonic oscillator ones with u — > 
w(f). 

\p\ + cu 2 (t)z 2 ] |nt >= (2n + l)w(t) |n* > (32) 

< ^ >= V^W = (^) 1/4 ^^H n (z^))e-^ 2 . (33) 

Since the creation and annihilation operators (&' (i) and in eq. fl2T?|) are constructed from 
2 and p 2 operators the eigen state \n t > of the number operator n(t)(—b' (t)b(t)) commutes 
with any function f(t) which is independent of z and p z . Therefore we find 

\p 2 z + oj 2 (t)z 2 + X j E(t)]\n t >= [(2n+l)w(t)+Aj£7(t)]K > . (34) 



By using the orthonormality relation of Hermite polynomials we find 

dz\Mz)\ 2 = Jdz\<n t \z>\ 2 = J dz{^)^j^\H n {z^{t))\ 2 e-^ 2 = 1 (35) 



dztf}^(z)ip m (z) — dz < n t \z >< z\m t > 



and 



^(^) 1 / 2 ^^(^^)i/ m (^V / ^))e^ ( ^ 2 = 5 nm . (36) 
By using / dz\z >< z\ — 1 in the above equation we find 

< n t \m t >= 5 nm . (37) 

The completeness condition is 

J2\n t ><n t \ = l. (38) 

n 

n our case the harmonic oscillator [— e 2 E 2 (t)z 2 +p 2 ,} which appears in eq. (127]) has imaginary 
3), Q| frequency ui(t) = ieE(t). Inserting complete set of harmonic oscillator states \n t > 
(J2 n \ n t >< n t \ = 1) in eq. (1271) we find 

—^22^2 J dt J dp J dp' Q J dp'o J dp'o J dp z J dq z J dz x J dz 2 J dz 

o it P0e iz Pz < po \ e -Ek£ P z\ p > Q > e -iz iPz < < ^| nt ><: ^isl-eSE^+pl-iXjEit)]^ > 

< m t \z 2 > |pg > e iZ2q * < p'o\e^^ qz \p^ > e^e"^" 

= l^J dt J dp J dp J dp J dp J dp z J dq z J dz 1 J dz 2 J dz 



1 d ■ 



-00 

e itp 0e izp z < p \ e -EkTtP*\ p > o > e -«U». < p > o \ < > e s[eE(t)(2n+l)~X 3 E(t)] < ^ > ^ > 

= 7^~j3 H / dt J dpo J dp' J dp'o J dp" J dp z J dq z J dz x J dz 2 J dz 

e it POe iz( Pz - qz ) < po \ e -^&P*\ p f o > e -iz lPz < p ^ < ^ > e sX d E(t) e -seE(t)(2n+l) < ^ > | p « > 

e i229z < Po |e^T | p ' " > e -«Po". (39) 
As advocated earlier the z integration must be performed from —00 to +00 for the shift 

n 

theorem to be applicable [5|]. Performing the z integration we get 

/ + °° dze iz ^~ qz) = 2rt5(p z - q z ) (40) 



and we find 

Fi(s) = 7A2E j + 2 dt I dp ° J dp °J dp °I dp ° J dpz J dzi J dz 2 e ltpo <po\e~^^\p' > 
e -iziP* < p / Q | < ^ > eS A^(t) e - se £(t)(2n+i) < > | p » > e i 22P2 < p »\ e -m&r*\p'>> > e-**&". 

(41) 

Inserting further complete sets of \t > states we obtain 

j^pE / dt J dtl J d h J d Po J dp'o J dpi J dpi' J dp z J dzi J dz 2 e %tpa 
< p \e~^^ Pz \p' Q > e~ iziPz < Pa\h >< h\ < z x \n t > e ^fi(t) e -.eE(t)(2n+i) < ^ > 
\t 2 >< t 2 \p'o > e iZ2P * < Pole^^bo" > e~ itp ° 

= H J dtjdtij dt 2 J dp J dp' Q J dpi J dp'o J dp z J dz x J dz 2 



-00 

1 d , 



e itpo < p \e~^)^ Pz \p' Q > e - iZlP *e- itlP '° < t x \ < z x \n t > e ^m e -seE(t)( 2n +i) < ^ > 
\t 2 > e lt2P °e iz * p * < Po|e^T^ Pz |p' " > e~ ltp " . (42) 

Since < n t \z > does not contain any time derivative ^ operator (see eq. (133]) ) we find (by 
using < ti\t 2 >= S(t x - t 2 )) 

Fj(s) = jt^YI J dt J dtl J d Po J d Po J d Po J dp'o J dp z J dzi J dz 2 

e itpo < p \e~^^ Pz \p' > e - izip *e- itip '° < Zl \n tl > e >*iWi) e -«*i(ti)(?n+i) < > jW> 

e iZ2P * < pi \ e m Pz \p% > e- itp ° . (43) 

1 d 1 d 

Since < po\e dt \p' > and < p'l\e E ^ Wo > are independent of the time derivative op- 
erator -| (they only depend on c-numbers p , p' Q , p'l and we can take < Po'|e'^w n * Pz |po" > 
and e~ ltp 'o to the left. We find 

Fj(s) = j^-^Yl J dt J dti J dp J dp J dpi J dpi' J dp z J dz x J dz 2 

e it( P0 -P'J>) < j%\e-m&'\p'» >< p \e~^^\p' > e - iziP *e- itlP '° < z x \n tl > e «W0 e -*^(*0(2»+i) 

< n tl \z 2 > e Ulp °e tZ2P \ (44) 

1 d 1 d 

It can also be shown that < p \e B w dt \p' > and < pl\e^ dt \p'l' > are independent of t. 
This can be shown as follows: 

< Po\f(t)j t \p' >= J dt' J dt" J dp'l" < Po \t' >< t'\f(t)\t" >< f |jC >< Po'\j t \Po > 
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= J dt' j dt" J dp™e- a 't>° 5(t' - f)f(i!')e*W ip' 6(j%" -p' ) = ip' I dt' e^ ^ /(f) 

(45) 

which is independent of t. Hence we can easily integrate over t in eq. ( 1441) . 

Integrating over t in eq. (jHJ) by using f^™ dte lt ( po ~ p °'^ = 2ir5(po — Pq"), we find 

Fj(s) = / dt ± j d Po J dp' J dpi J dp z J dzi J dz 2 < pg|e^^ p *|po >< Po\e~^ ipz \p' > 

e -izi P z e -it iP > < > eS A J B(t 1 ) e -se J B(i 1 )(2n+l) < > gi^'gfaap, _ ( 4 g) 

Since J dpo|po >< Pol = 1 and / dp z e iPz ( Zl ~ Z2 ^ = 2n5(zi — z 2 ) we find 

Fj{s) = / d Pto/ dz iJ dz 2 e~ lzlPz < Zl \n tl > e »^(tO e -«E(tO(2n+i) 

< n tl |z 2 > e iZ2P2 = -^-V / dt / dp' / dz\ < z\n t > \ 2 e -A j E(t) e -8eB(t)(2n+i)_ ( 47 ) 
(27r) n J J J 

Since the harmonic oscillator wave function is normalized (see eq. ( 1351) ): 



J dz\ < z\n t > | 2 = 1 (48) 



we find 



Using the Lorentz force equation 

dp^ = eF^ u dx u (50) 
we find (when E(t) is along z-axis, eq. f|T4l) ) 

dp = eE(t)dz. (51) 

Hence 

FJs) = -±- dt dz \ (52) 
JV ; (4tt) 7 J smh(seE(t)) v ; 

Using the above expression for F in eq. ( 124"1) the effective action becomes 
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Using the values of Xj (the eigen values of the Dirac matrix 7°7 3 from eq. (fl9l) ) we find 

= J_ f°° ^£ / fa f d 2 p T e is{p T +m2+ie \eE(t) coth(seE(t)) - -]. (54) 
871- 3 Jo s J J s 

The imaginary part of the above effective action gives real particle pair production. Using 

the series expansion 

+ (55) 



sinhx x "^7r 2 n 2 + x 2 

at a 



the s-contour integration is straight forward [1|, [3|, [6j. Performing the s-contour integration 
around the pole s = r^ji we find for the probability W of pair production (which is twice 
the imaginary part of the effective action) 



W = 2ImS (1) = — J] - / (fx / d 2 p T |e^(t)|e" n7r ^W. (56) 
n= i w ■/ ■/ 

Therefore the probability of producing a fermion (or antifermion) per unit volume per unit 
time with transverse momentum pt from arbitrary time dependent electric field E(t) is given 
by 

dW \eE(t)\ T f1 y T +m2 , 

- ! Log[l - e ^i^mr] (57) 



d 4 xd 2 pT 47T 3 
which reproduces eq. (jSJ). 

To conclude we have obtained an exact one-loop non-perturbative result for the probabil- 
ity of fermion-antifermion pair production per unit time per unit volume per unit transverse 
momentum (of the fermion or antifermion) from the arbitrary time dependent electric field 
E(t) via Schwinger mechanism. We have found that the exact non-perturbative result is 
independent of all the time derivatives d f„ , where n = 1,2, ....00. 

The quantum back reaction problem for the Electric Field is studied in [8|. There it is 
shown that the exact numerical solution is consistent with the solution of the semi-classical 
transport equation having a source term for pair production which utilizes Schwinger's 
constant electric field E result but with the replacement E —> E(t). This result seems to be 
consistent with our finding here. However in those simulations the Schwinger source term 
was modified by Bose Enhancement effects for boson pair production and Pauli-blocking 
effects for fermion pair production 8|, |9j. These Pauli-blocking effects are not present in our 
calculation which needs to be better understood. Particle production from arbitrary time 
(and space) dependent background fields is important in early universe dynamics and in the 



production and equilibration of the quark-gluon plasma at RHIC and LHC 
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